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. . . We display an interesting sum rule for the dynamical thermal conductivity for many standard 

I models of condensed matter in terms of the expectation of a thermal operator. We present the 

. thermal operator for several model systems of current interest, which enable an evaluation of the 

' sum rule and the Lorentz number, the thermo electric figure of merit as well as the thermopower at 

high frequencies. As a by product, we present exact formulae for the T = chemical potential /i(0) 
' for charged many-body systems, including the Hubbard model, in terms of expectation values of 

^ ' extensive operators. Simple estimates are provided for the thermopower of an infinitely correlated 

band model on the triangular lattice, modeling the physics of the sodium cobalt oxide system. The 
' present result goes beyond the Heikes Mott formula for the thermopower, and contains an additional 

transport correction that is sensitive to the lattice topology as well as the sign of hopping. 
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I. INTRODUCTION 



A problem of considerable current interest is the evaluation of the thermal response functions; k the thermal 
' O ■ conductivity and S the thermoelectric response function (Seebeck coefficient) for complex systems. These include 
g strongly correlated matter, where the role of interactions is profound and impossible to capture within perturbative 
^ formulations, and also materials with complex unit cells, such as the skutterudites [1], and negative thermal expansion 
systems where geometrical frustration plays a major role [2]. One major motivation for these efforts is from the 
requirements of the thermoelectric industry, where a large thermoelectric figure of merit (defined in Eq(39)) is desirable 
^ , , and it is important to understand the material characteristics that promote a large figure of merit. The standard 
y—i • tool, namely the Boltzmann equation approach, is overly restrictive in that it is tied down to the concept of long lived 
^-H ^ and weakly interacting quasi particles, which are less useful in these contexts. 

In this paper, we present a new formalism for computation of the thermal response functions, by considering 
, the response to dynamical temperature gradients, i.e. a frequency dependent problem. This yields the frequency 
' dependent thermal conductivity k,{lij) as well as the thermopower S{lo). The frequency dependence of k,{lj) is of 
] relevance theoretically and also practically: a silicon chip operating at several gigahertz clock speed in a computing 
■ processor needs to transport the AC Joule dissipation. The relevant transport parameter in this context is the 
] dynamical k,{uj), which can be much smaller than the static k [3] for uj ^ lOOGHz. 
^ ■ Although studying the dynamical conductivity looks at first sight like further complicating an already intractable 
' , . problem, in fact it turns out that several simplifications arise in the limit of high frequencies. These have a close 
^ ' counterpart in our fairly successful theory of the Hall constant of strongly correlated matter [4] . In brief, the problem of 
Q ' transport may be decomposed into the following conceptually distinct but technically intertwined aspects (a) transport 
O ■ proper (b) single particle description and (c) interactions. The single particle aspect determines such variables as the 
^ ] density of states and its derivatives, the transport aspect involves the elastic and inelastic scattering rates usually 
modeled by a relaxation time, and the interactions aspect involves the influence of correlations. While none of these 
. , are trivial, the last one, namely interactions is conceptually as well as computationally most difhcult. For example 
' in Mott Hubbard systems, the very notion of the hole is non trivial in the context of say the Hall constant, where 
one knows that near half filling, the momentum space definition can give opposite results in certain cases to the real 
space picture for the sign of the Hall constant [4]. 

In this context, going to high frequencies is a great advantage if one can isolate combinations of thermal response 
functions that are relatively less frequency dependent. In the Hall constant problem, the Hall resistivity p^^ is quite 
benign, being w independent at least in the Drude theory, as opposed to the Hall conductivities a^^ and tr^^, which are 
serious functions of w individually but not in the combination leading to p^^ . Moreover the high frequency limit is the 
starting point of the Mori type treatment of transport quantities, as seen in the example for the Hall constant [5] . In 
a similar spirit, we identify several variables in this work, and write down the operators that need to be evaluated. We 
stress that the expectation values of the objects computed here are equilibrium values, and while quite non trivial in 
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detail, arc conceptually much easier than the transport objects- this is the advantage of disentangling dynamics from 
interactions as in other contexts. While our approach is not exact either, as it depends on the assumption of weak 
frequency dependence of such combinations, it is quite distinct from the Boltzmann approach and provides a useful 
counterpoint to the latter, treating the effect of interactions more respectfully. In simple examples, we show that the 
weak Lj dependence is an excellent assumption, and further our results should provide a stimulus for experimental 
investigation of the frequency dependence of thermal transport constants. 

We present an interesting sum rule for the real part of the dynamical thermal conductivity. The sum rule is in 
terms of the expectation of an operator that is presented for several models of condensed matter physics, and while 
the task of evaluating it is non trivial in detail, it is quite feasible to estimate it from other considerations that we 
present. The situation is akin to the well known plasma or the f-sum rule; the Thomas Ritchie Kuhn sum rule gives 
the integral of the real conductivity in terms of easily measurable variables in the continuum Airnq^/m, whereas on a 
lattice, it is (see Eq(All)) the expectation of essentially the kinetic energy [6]. Such an expectation is non trivial, but 
several numerical techniques lend themselves to this task quite effectively, such as exact diagonalization and Quantum 
Monte Carlo methods. This new sum rule should be of considerable interest in many situations, including the case 
of non linear lattices, where the classical counterpart is already interesting, since it is the expectation of a single 
extensive operator, and apparently equally unknown. 

Our paper has a twofold purpose, we present the formalism for the thermal transport variables in some detail with 
the hope of returning to calculate some of them in a later work, and the other is somewhat pedagogical in nature. 
Existing literature is quite focusscd on specific issues, and many subtleties receive less mention than their due. It is my 
impression that a detailed description of the formalism might be useful to some workers. The results of computations 
based on the formalism presented here should shed light on currently interesting issues, such as the departure from 
Wiedemann- Franz ratio for the case of strongly correlated systems as seen in experiments, where the Lorentz number 
departs from the free electron value at low T [7]. 

An interesting by-product of our work is the exact formula for the chemical potential /i(T) at T = for charged 
many body systems in terms of the expectation of an extensive operator. In general the only way of computing /i(T), 
including in the ground state, is from the thermodynamic formula /x(T) = ^'^j^^^ ■ However at T = 0, our result 
shows an alternate and potentially widely useful method. This arises from the physical requirement of the vanishing 
of thermopowcr at all frequencies in the groimd state, being related to entropy transport, and since the high frequency 
formulae explicitly involve the chemical potential, these yield the new formulae for the zero temperature chemical 
potential of all charged many body systems. We comment on this further in the discussion section (X) . The plan of 
the paper is as follows. 

• In Sections II and III, we derive the complete expressions for the finite frequency thermal conductivity and ther- 
mopower, using standard linear response theory, and obtain a sum rule for the dynamical thermal conductivity. 
It is expressed in terms of a new extensive object ©^^ that we term as the "thermal operator" . 

• Section IV makes a connection with the "standard" Kubo formulas. Our results contain corrections to the Kubo 
formulas for the case of non-dissipative systems such as superconductors, and it might be of interest to see the 
origin of these in the fine tuning of the standard derivations, as provided in Section IV. 

• In Section V, we identify the high frequency observables that are suitable for study being less sensitive to uj. 
The high frequency Seebeck coefficient is amongst these, and is expressed in terms of another new extensive 
object that we denote as the "thermoelectric operator" . 

• The thermal and thermoelectric operators have a dependence on the interactions as well the details of the specific 

model systems. Hence it is important to compute these operators for standard models, which is accomplished 
in the latter sections. In Section VI, we present the thermal operator Q^^ for the anharmonic disordered lattice 
vibration problem, and make a connection with the classical limit. Illustrative examples from harmonic lattice 
are given to provide a feel for the nature of the operators encountered. 

• In Sec VII, we study the intermediate coupling models of current interest, namely the Hubbard model, the 
homogeneous electron gas and the periodic Anderson lattice, for which we present the thermal and thermoelec- 
tric operators, and also indicate the role of disorder in modifying these variables. We establish the essential 
correctness of this approach for the case of zero coupling, where the standard Boltzmann- Drude results for the 
Lorentz number and the thermopowcr are reproduced. 

• In Sec VIII, we study the thermal and thermoelectric operators for the strong coupling models: the Heisenberg 
model and the infinite U Hubbard model. 
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• In Sec IX we apply the results of Sec VIII to the currently interesting case of the triangular lattice NaxCo02, 
a metallic thermoelectric material. Our theory yields a formula for the High T limit of the Seebeck coefficient 
for this system, which contains the Heikes Mott formula and provides a transport correction to the same. The 
large transport correction highlights the role of the lattice topology. 

• In Sec X, we discuss the special nature of the thermal type operators, namely the vanishing of its expectation 
in the ground state. We also comment on the exact formulae for the chemical potential at T = for charged 
fermi systems that our approach yields. 

• Appendix A contains the already well known results for the electrical conductivity for completeness. 



II. THERMAL CONDUCTIVITY AT FINITE FREQUENCIES 

For models involving particle flow, we recall that the heat current is defined as the energy current minus times 
the particle current. Therefore to generate the linear response equations, we write the grand canonical ensemble 
Hamiltonian in the presence of a temperature variation as (wc = w + iO+) 

K = + K^e-^'\ (1) 

with adiabatic switching from the infinitely remote past t = —oo as usual, and Kq = ^fK{r) = ^f{H{f) — ji n{r)). 
Here H{r) is the energy density, and since we are mainly dealing with lattice models, we sum over r. The operator 

K^ = Y,mK{^, (2) 

r 

where ip{r) is a small (pseudo) gravitational field with some spatial variation such that its average is zero. This 
expression follows the conceptually important work of Luttinger [8] . Alternate schemes such as the Kadanoff Martin 
method [9] lead to the same results. Here ^/;(r) is not quite a gravitational field since it couples to K{r) rather 
than to -ff(r), but serves the same purpose, it is a mechanical field as opposed to the true spatially dependent 
temperatiirc, which is a thermodynamic field. The latter is more subtle, and ensuring that such a varying field is 
established in a given system is in general a very hard thing to achieve rigorously- requiring as it does a detailed 
understanding of the equilibration processes that operate on a microscopic and inhomogeneous level. Luttinger's idea 
reorganizes and subdivides this complication into logically separate parts, and allows us to deal with a mechanical 
field with an arbitrary spatial variation. The connection with true temperature variation is the next logical step, 
where Luttinger showed that we may essentially regard V'!/^(r^ = ^^r^. Here we denote the local temperature by 
T(r), and /3(r) = ^^^^^ A formal discussion is provided by Luttinger [8], who drew an analogy with the Einstein 
relation. The latter relates the response to the (mechanical) electrical field to the observed concentration gradient, 
which in turn is related to its conjugate namely the (thermodynamic) chemical potential gradient. The resulting 
frequency dependent linear response functions are understood as responses to this mechanical force. Relating these to 
the thermodynamic response functions is the other part of the issue, and may be proved only under certain standard 
assumptions [8] ; we are content to assume the same in this work. A simple qualitative argument to help motivate the 
identification Vip{r) = is that the temperature profile varies in such a way as to annul the variation of the 

added energy fluctuation, therefore if locally tp{r) increases so does T(r) so as to maintain the local constancy of the 
action J2rf^i^)^i^- 

Let C be the number of sites, or equivalently the volume of the system (by setting the lattice constant as unity). 
We next fourier transform Eq(2) by introducing 

V'(r) = '4>{k) exp {—ik.f) (3) 

fc 

= 7 E ^(^) exp(-tfc.f), (4) 
k 

so that 

K^=Y,K{-k)m- (5) 
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The heat current operator is Fourier decomposed as J'^{r) = ^ exp (— zfc.r) J'3(fc). The induced heat current in 
response to the fcth mode with can be computed from hnear response theory. We will assume a "temperature" gradient 
along the unit vector x, so that VT(r) = TVtp{f} = T^{—ikxx)'4>{kx) exp {—ikxx). Then 



K,{u)c) = lim ^(fca;, o^c) with 

fea;— ►o 



(6) 



The expectation {j§{kx)) is given by standard linear response theory [10,8,11] and using Eq(6) we find 



K{kx,OJc) 



= Wk^l /I dt'{[jS{kx,t),k{-kx,t')]) 



fuVcTkxjC 



{[J^{kx),ki-kx)])+i / e-=*' dt'{[J^ikx,t'),[k{-kx,0),K]]) 



(7) 
(8) 



We integrated by parts to get the second line from the first. In the limit of almost uniform variation, kx — > 0, we 
show that 



{[jQ{kx),k{-kx)])-^-kx{en, 



(9) 



where is the thermal operator, to coin a name. In the uniform (i.e. kx — > 0) operator is not a constant of 
motion, and yet the average is of 0{kx), this is so since the uniform term vanishes by noting that the thermal average 



(Identity I), 



with Z = Tre , the last identity following from the cyclicity of trace. Therefore we can write an equation for the 
thermal operator directly as 



e'^'^ = -ihn^—[jS{kx),k{-kx 



(10) 



and thus it is straightforward if tedious to compute it, given the explicit forms of the current and energy operators. 
The heat current is obtained from the continuity equation for heat density, written in momentum space as 



\im^j-[k{-kx),K]=nj^, 
where = J§{k)/k^o, so that the thermal conductivity at finite frequencies is 



k{lUc) 



1 



JkOrT C 



POO 

/ e'^^*' dt'<[jQ{t'),jQm> 
Jo 



(11) 



(12) 



It is worth making a comment on the difference between this calculation and that of cr(w), the electrical conductivity 
at this point. A calculation of a{oj), as in Appendix A, is in close parallel to this one for k{oj), with the electrical 
current and charge density fluctuation operators (deflned in Appendix A) replacing the heat current and the energy 
density as in J^{kx) Jx{kx) and with K{—kx) p{—kx). In that case, the commutator [Jx{kx)i p{—kx)\ explicitly 
vanishes as kx 0, and hence the analog of the first term inside the round bracket of Eq(8) automatically begins as 
0{kx), with a coefficient oc (r^^) ( defined in Eqn(A4)). In the present case, Eq(9) together with Eq(ll) again leads 
to a cancellation of kx between the numerator and the denominator, leading to a finite result in the uniform limit. 
Thus the thermal conductivity calculation is somewhat disguised by the non-commutation of the uniform operators 
i.e. [J^, K] 7^ 0, which is luckily immaterial, since the expectation of this object vanishes due to the Identity I. 
We perform a Lehmann representation [11] to write Eq(12) with p„ = exp(— /3e„) as 



k(Wc) 



i 1 



Pn -Pn 



Cm + hAJ, 



■\{n\jQ\m)\' 



(13) 
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Using the partial fractions identity (for any A) 

1 



1 / 1 



we obtain 



k{u1c) 



hu)c{hu)c + A) A \TiuJc hojc + Ay ' 



(14) 



where 



Dq = ^ 



Pn -Pn 



■\{n\J^\m)t 



1 

£ 



Jo 



(15) 



Using the standard identity l/(wc — E) = P — E) — in5{w — E) for any E, it is instructive to write the real 
part of the Lehmann representation of Eq(14) as 



Re k{uj) = — (5(a;)DQ + Re Kreg{i^) with 



Re K,reg{oj) = — ( ) 2_, Pn\{n\J^\m)\^ 6{em -en- hul), 



TC 
1 



.|(n|J«|m)|^ 



(16) 
(17) 



Here Kreg is the regular part of thermal conductivity that excludes the delta function at w = 0, and Dq is the weight 
of the delta function at w = 0. It is seen that the weight of the zero frequency delta function involves a different 
object than the one in Eq(15), there is a cancellation of terms with equal energy between the two terms in Eq(13). 
We elaborate on this distinction below. 

The sum rule for the real part of the thermal conductivity (an even function of w) follows from Eq(16) by integration, 
or cquivalcntly from Eq(13) and using the Kramers Kronig relation to identify the coefficient of i/to at high frequencies 
in the complex conductivity as the integral of the real part: 



f 

Jo 



Re K{u})du) = 



TT 



2hTC 



(18) 



This interesting sum rule is one of the main formal results of this paper. The extensive thermal operator ©^^ is 
analogous to the stress tensor r^^ in the sum rule Eq(All) for the electrical conductivity a{oj), and depends in its 
details upon the underlying model; we present it in the case of several problems of current interest in the following. 
We may rewrite Eq(14) in a more compact form as 



(19) 



The first term in Eq(19) is oc Dq, and is in addition to the original Kubo formula [10] for thermal conductivity, namely 
the second term. It represents non trivial a correction for non-dissipative situations, but vanishes in dissipative cases 
due to reasons that we elaborate below. 

Integrating the second term of Eq(19) gives an estimate for Re K{u))duJ namely ^^^^ dT{J^{—iT)J^{0)), 
which has been written in literature earlier in its high temperature limit [12]. However it cannot be viewed as a sum 
rule, since the estimate is a correlation function of a pair of extensive currents, rather than a direct expectation value 
of a single extensive operator as in Eq(18). Further, in non-dissipative situations, this estimate misses the contribution 
from the Dq term and is thus incorrect. For dissipational systems, it is does reduce to Eq(18), on using the vanishing 
of Dq and Eq(15). Therefore it is clear that in all cases, Eq(18) stands in complete parallel to the f-sum rule Eq(All) 
for (t{u)) with the thermal operator ©^^ relacing the stress tensor r^^. 
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Wc next make a few comments on the object Dq, which enters into the correct expression for thermal conductivity 
Eq(19), and Dq that enters in Eq(17). Note that the difference between Dq and Dq lies in the nature of the double 
sum over the current matrix elements, Dq explicitly excludes terms with equal energy, so that the difference is a sum 
over all degenerate manifolds (including diagonal matrix elements n = m): 

Dq-Dq = -np J2 Pn\{n\J^\m)f (20) 

A detailed discussion of the charge stiffness in the similar context of charge transport is given in Ref [13], where the 
distinction between zero and finite temperatures is elaborated upon. Our discussion here is along similar lines, except 
that here we use the term charge stiffness to denote both Dc and Dc whereas these are denoted by the terms Meissner 
stiffness and charge stiffness in Ref [13]. In the case of the thermal conductivity Dq and Dq are thus the thermal 
analogs of the Meissner stiffness and the charge stiffness of Ref [13]. 

The value of Dq is zero for most dissipational situations, i.e. for generic systems at finite temperatures. It is 
however non zero in general for systems that display macroscopic many body coherence, i.e. either superfluidity or 
superconductivity. For a superconductor, the analogous object, D^ in the charge conductivity (see Eq(A7)) is in fact 
related to the London penetration depth or the superfluid stiffness, and its non vanishing is the very hallmark of its 
"super" -nature. It is alternately obtainable from the Bycrs-Yang relation D,, cx lim^^o ^^-^(f/*) where F{<j)) is the 
free energy of the superconductor in the presence of a hole threaded by a flux </>. While there is no obvious parallel 
to the flux for thermal conductivity, the formal expressions are very similar. 

There has been considerable discussion of the corresponding charge stiffness in recent literature [14] for integrable 
many body systems in low dimensions, such as the Heisenberg or the Hubbard models in 1-d, where the existence 
of many conservation laws pushes them away from equilibrium with negligible restoring forces, leading to a ballistic 
rather than diffusive behaviour, as first pointed out by Giamarchi [14]. An interesting numerical study of the statistics 
of the Kubo conductivity in 1-d has also given some insights around the zero frequency limit [15]. Very recently the 
thermal conductivity has also been discussed in literature for integrable models [16]. In general terms, integrable 
models are akin to the free particle system rather than to superconductors or superfluids. since there is no broken 
symmetry. In such cases, we expect Dq = but it is possible that Dq ^ 0, albeit possibly only for finite sized 
systems. This is analogous to the charge stiffness story [13], where D^ is zero for normal metals whereas Dc may be 
non zero. In this sense the zero temperature limit is singular; Dc occurs in both conductivity and in the Byers-Yang- 
Kohn type formula Dc cx lim^^o 's^^o{4>)i with EQ{(j)) the groimd state energy in the presence of a flux [13], whereas 
at non zero T there is a distinction as in the RHS of Eq(20, A13). The object Dq involves matrix elements of the 
current in states with distinct energy, and these arc zero if the energy/heat current is conserved, as in the case of 
the free electron system. More non trivially, in the case of the exactly integrable Heisenberg model in 1-d, the higher 
conservation laws include the energy current [16]. In such a case, since the heat current is a constant of motion, it 
follows that J^{t) = J^(0), leading to a simple calculation of K(a;) using Eq(12). Here the second term in Eq(12) 
vanishes on using J§{t) = J^(0), giving 

"("'J = sir*''"'- 

Thus k{uJc) may be computed for integrable models with a conserved energy current directly from a knowledge of 

We will show in Sec IV that the vanishing of Dq is natural for generic systems, using an argument based on the 
Kubo identity, and also pinpoint the possible technical reason as to why it is non zero for say a clean superconductor. 
In subsequent sections we provide an evaluation of the thermal operator O^^ given in Eq(lO) for various models of 
current interest. 



III. THERMOPOWER 

We next study the linear response formulation for thermopower. If we define the charge current and its Fourier 
decomposition as J(f) = ^ ^exp(— ifc.r) J(fc), then in the presence of an electric field and a temperature gradient: 

< >= a{iv)E^ + 7(a;)(-VT), (21) 

where Jx = Jx{kx — > 0). The thermopower is defined as 
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5H = 



cr(a))" 



We can obtain 7(0;) by linear response theory in parallel to the thermal conductivity, and the answer is 



7(wc) = lini 

The expression for {Jx) follows the same lines as that for the heat current and we find 



Jo 



where we have introduced another important extensive operator, the "thermoelectric operator": 



(22) 



(23) 



(24) 



It is useful to perform its Lehmann representation, so that 



7(wc) = 



TVjJrTC 



Pn -Pr, 



\l\m){n\J^H 



(25) 



(26) 



We can proceed further as in the case of thermal conductivity by expanding the product using partial fractions and 
find 



7(wc) 



Pn-Pm \ {n\X\m){n\J^\m) 

huJcT ' ' TC ^ \ Cm - en J Cn- em + TVj^c 



^^7 + ;^E 



(27) 



where 



Vn Pm 



\X\m){n\J^\m) 



(28) 



Since the two currents involved are not identical, it is not possible to take the real part of this expression readily, 
and it is not profitable to seek a sum rule. It is however, interesting to note that by combining various formulae, the 
thermopower has a high frequency expansion of the type 

S(uj) =S* + 0(-) where S* = (29) 

and so a knowledge of ^^'^ is useful in determining the asymptotic behaviour. We can rewrite this equations as 



7 1 Z""" - 



7(wc) 



We again note that = for most dissipational cases. 



IV. KUBO TYPE FORMULA FOR GENERAL NON DISSIPATIVE SYSTEMS 



(30) 



In an earlier section, we have derived the linear response theory expressions for the various conductivities Eqs(6), 
and argued that the usual formulae quoted in literature e.g. in Ref [17] following Kubo [10] and others [8] have extra 
terms that arise for non dissipative systems, such as superconductors and superfluids. The reader might wonder as 
to the origin of these extra terms, since the common belief is that the Kubo formulas are complete and exact. In 
this section we obtain these terms from a careful analysis of the Kubo arguments. We deal with the case of thermal 
conductivity below, but the arguments apply equally to the electrical conductivity as the reader can easily see. In the 
case of the electrical conductivity a simpler argument gives the final result using the transverse gauge, i.e. by realizing 
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the electrical field as the time derivative of a vector potential rather than as the gradient of a scalar potential [18,19], 
although the present method is readily generalized to that case as well, as shown in the Appendix. 
We begin with Eq(7) rearranged in the form preferred by Kubo ; 

^{k^cOJc) = ^^j^ dt e"^^'Tr {[po,K{-k,,-t)]J^{k,)) , (31) 



-I3K 



with the density matrix po = ^ ^ • From this point onwards one uses three important steps 
1. The so called "Kubo identity" [10] which is written as )C[K{—kx, —t)] = where 

r0 



)C[A] = [po,A]- f po[A{-iT),K]dT, (32) 

^0 



which is easily derived by inserting a complete set of states. 
2. The conservation law 



^[k{-k,,-t),K] = jQ{-k,,-t). (33) 



3. The homogenous or uniform limit k^ 0. 

These three steps together give the usual stated result Eq(12) without the first term involving Dq. These steps lead 
one to an expression involving the matrix elements of K{—kx, —t)/kx and it is tempting to replace this with those 
of J^(0, —t). It is precisely here that one needs to be careful. As long as kx is finite it is clear that K{—kx, —t) as 
well as j§{—kx, —t) have matrix elements between eigenstates of the Hamiltonian K that are at different energies, 
i.e. inhomogeneous excitations cost energy. In the limit of fcx — > this need not be so, the current operator can have 
matrix elements between states of the same energy, whereas by construction the commutator [K{'-kx, ~t), K] filters 
out states with the same energy. Thus the operator identity Eq(33) is potentially in trouble precisely where we need 
it. What is easier to justify is to use a time derivative to filter out equal energy states so that a good alternative 
identity to Kubo's identity is to use 

dlC[K{-kx,-t)] ^ ^ 

dt ' ^ ' 

This is the Kubo identity again, but with the operator (— i)[_R'(— A-j;, — f), if] in place of K{—kx,—t), this extra 
commutator with K is transferred to J^{—kx, —t) on using the conservation law, and improves matters. In this form 
we can take the next two steps safely since [J^(0, —t),K] does not have any matrix elements between equal energy 
states. We can now integrate the Eq(34) between finite times and get 

K[K{-kx, -t)] = K[K{-kx, 0)], (35) 

This equation tells us that the errors in replacing the matrix elements of K{—kx, —t)/kx with those of J^{—t) are 

time independent, i.e. relating solely to the matrix elements of current in the manifold of zero energy difference states. 
Such a manifold of states is statistically significant only in a superconductor, and hence we expect it to be relevant 
in that context. 

Proceeding as before to the uniform limit, we get 

1 /■ oo P0 

^Tcl "'^^''^^ Jo '^^(JS(-t-i^)J^(^))- (36) 

The term in curly brackets is recognized as CDq as in Eq(19), and hence the two derivations lead to exactly the same 
answer. 

We have seen that Dq owes its possible non vanishing to the matrix elements of the current operator in the manifold 
of zero energy difference states, such as exist in a superconductor. For most dissipative systems, such a manifold is 
statistically insignificant, and it should be noted that in such instances, Dq is actually zero. 
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The only instances where Dq and analogously Dc of Eq(15, A7) are non zero have to do with supcrfluid or 
superconducting states. We remind the reader of the change in terminology from Ref [13] as remarked on earlier 
below Eq(20), Dc represents the Meissner stiffness of Ref [13]. In these cases, Dc for example is proportional to the 
superconducting density of electrons as in the London theory. Traditionally ps is obtained by taking the current 
current correlator and taking the limit fc — > for the transverse part, i.e. J.k = 0, after having taken the static limit 
at the outset- this is the classic Meissner effect calculation [20] . However, it is also natural to see the superconducting 
fraction in the conductivity where wc first set fc ^ and then let lo ^ 0. This is the procedure followed in the 
well known Ferrell-Glover-Tinkham sum rule determination of the penetration depth in superconductors from the 
infrared conductivity [21]. In this sequence of taking the uniform limit, it should not matter if we have worked with 
the longitudinal conductivity rather than the transverse one, and our derivation here is the longitudinal counterpart 
of the transverse derivation first obtained for the Hubbard type models in Ref( [19]). A nice discussion of the various 
limits is given in Ref( [22]). A superconductor would also have a non zero Dq if impurities are neglected, however 
impurities do scatter quasi particles and hence one expects disorder to make this object tend to vanish, unlike the 
variable Dc which is less sensitive to disorder by virtue of the Anderson theorem. There is also the similarity of the 
Dq term with that of second sound that one expects in a superconductor, namely an undamped energy wave; Dq 
then appears as the residue at a second sound pole in the uniform limit of a suitable propagator. 

This derivation also indirectly shows that in the normal dissipative cases, when Dq = 0, the sum rule Eq(18) is 
obtained starting from the conventional Kubo formula, i.e. Eq(19) without the first term, by integrating the current 
current correlator over frequency. This is so since such an integral gives the second part of the expression for Dq 
which must equal the first for it to vanish! 

It is worth mentioning that another independent argument for the vanishing of Dq and likewise Dc, or in another 
context, the spin wave stiffness of a fully saturated ferromagnct [23], uses the notion of a position operator. If a 
position operator, or a similar one for energy transport such as the boost operator ( J drfH(r)) is allowed in the list 
of allowed operators, then it is easy to show that Dc type operators are zero, since one can rewrite the current matrix 
elements in terms of the matrix elements of the position operator, and then the second part of Eq(A7) becomes the 
commutator of the current with the position operator and hence cancels exactly the r^^ term. Such a procedure is 
invalid in general, as explained carefully by Herring and Thouless [23], the position operator is illegal in the Hilbert 
space where periodic boundary conditions are used in a continuum field theory. For a lattice field theory of the sort 
we are considering, we can always force the introduction of a lattice position operator X = '^^XjUp^, but now a 
careful calculation shows that the cancellation is incomplete, while the bulk terms do cancel, the boundary terms do 
not. 

From this line of thought we conclude that Dc and its relatives would vanish, if the nature of the quantum states is 
such that the error made in introducing the position operator is negligible, this again is plausible in a highly dissipative 

system, where phase coherence is lost over some microscopic length scale that is shorter than the system length. For 
a superconductor such is certainly not the case, and nor for that matter, in a ferromagnet Ref( [23]). 



V. HIGH FREQUENCY RESPONSE: LORENTZ RATIO AND THERMOELECTRIC FIGURE OF MERIT 

Of great interest are the two quantities, the Lorentz ratio 

and the thermoelectric figure of merit Z times T (a dimensionless number) 

Zi.)T='-^M^ (38) 



k{uj) 

We can readily evaluate these in the high frequency limit as 



(39) 



(0xx^ (.y-xa:^ 



l(hxx\2 

Z*T= , , / , . (41) 
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Let us also note that the thermopower at high frequencies is expressible in terms of the operators we displayed 
earlier; 



(42) 



These variables are computable much more readily than their dc counterparts, and since these do capture the effect 
of interactions, they are of considerable interest. 



VI. LATTICE THERMAL CONDUCTIVITY 



We consider here the sum rule as applied to the case of the lattice vibrations. This is a field of great current activity 
[24,25], where most studies are in low dimensions and treat classical anharmonic and disordered lattices with a view 
to study the conditions for existence of a Fourier law, and also for the emergence of a finite conductivity. We consider 
a general anharmonic disordered lattice with a Hamiltonian 



A. 

2m i 



(43) 



where Vij = Vj^i = V(ui — uj) is a symmetric two body potential that is an arbitrary function of the displacement 
variables Ui including possible anharmonicity. Since no particle flow is involved, it suffices to consider the canonical 
ensemble Hamiltonian. We will denote the equilibrium lattice points by Rj = {Xj.Yj, ..) and the antisymmetric force 

vector Fi_j = —j^Vij. In the following, wc will need the derivatives of the antisymmetric force vector which will be 



denoted by F. 



-F^-, and have the obvious symmetry property F', 



x{n) 



(— 1)" ^F^^^^\ The energy current 



can be obtained from computing the time evolution of the energy density, and combining it with the conservation law 
thereby giving 



(44) 



.-here E{k) = ^, 



Hj. This form of the current operator is popular in literature, e.g. Ref( [26]) and differs from 



the one given by Hardy [27] in some details, which are not expected to be very significant. To aid the passage to 
a classical limit we display the dependence on Planck's constant in this section, setting it to unity elsewhere. The 
current at finite wave vectors can be found from the above as 



{X, - X, ) 



rrii 



(45) 



where we have chosen k along the x axis. Using the various definitions, we find the thermal operator 



-I 2 



1 \ - (Xa - Xj)'^ 



(46) 



This together with Eq(18) then gives us the sum rule for thermal conductivity. Note that the classical limit exists 
and is readily found by dropping terms involving h in this expression. In the classical limit of Kubo's formula Eq(19), 
the inner integral over r collapses to give 



l^classicali'^c) — ^ J 6*'^"* dt {J^ {—i)J^ ifi)) classical i 



(47) 
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and it is clear that we can express the sum rule in terms of the expectation of {J^^x)- However our formula Eqn(18) 
provides a slight advantage in this case too: our result implies that 



n 



classical 



classical ; 



(48) 



a result that is non trivial since the LHS involves a single extensive operator whereas the RHS involves the square of 
an extensive operator, and moreover the only way to prove it seems to be to go back to the inhomogeneous (i.e. k 
dependent) energy density and perform the steps stated above! 

As a simple illustration consider the case 1-d Harmonic lattice, H = "^j^Pi / {^m) + ks/2'}2^{ui — u^+i)^, with 
ks = mwQ. In this case, we expect heat transport to be ballistic rather than diffusive, and hence the static thermal 
conductivity to diverge. Unlike the case of free electrons (discussed elsewhere in this paper), the current operator is 
not completely diagonal in the normal mode (phonon) operators, in contrast to the Peierls form of the energy current 
[28], namely = ^^^kVknk- The integrated conductivity is non trivial, it gets its weight at finite frequencies from 
phonon pair creation and destruction processes. The force is linear in the displacements and hence we work out the 
thermal operator: 



(49) 



The expectation of this object can be computed at finite temperatures on using the phonon harmonic oscillator 

representation as 



1 1 

^ 2 ef^'^" - 1 



)j'±cos{k) + ^sm\k) 



(50) 



where LUk = 2a;o sin(fcao/2) is the acoustic phonon energy. At zero temperature, i.e. in the ground state it can be 
shown to vanish exactly; this precise cancellation is a characteristic feature of the thermal operators that is in common 
with most other quantum systems, and we discuss its connection with the vanishing of the specific heat in another 
section (see discussion near Eq(89)). Similar expressions can be worked out in any dimension d. At low temperatures, 
it is easy to see that O^^ oc 7*^+^ in d dimensions, so that the sum rule Eq(18) is oc T"^, i.e. similar to the lattice 
specific heat. At high temperatures, since this model has no scattering of phonons, the sum rule increases linearly 
with T indefinitely. Through Eq(18) this implies a T independent behaviour at sufficiently high T. The computation 
for the non linear lattice would contain the effects of umklapp as well as normal scattering, and should be numerically 
feasible as well as interesting. 



VII. INTERMEDIATE COUPLING MODELS: 
FINITE U HUBBARD MODEL, HOMOGENEOUS ELECTRON GAS AND PERIODIC ANDERSON 

LATTICE 

In this section we will consider the popular intermediate coupling models, namely the Hubbard model and the 
electron gas as well as the Periodic Anderson Lattice. We derive for these models the Thermal operator G^^and 
the Thermoelectric operator We further demonstrate the formalism in the limit of zero interaction, where 

the Boltzmann Drude results are reproduced for the specific transport variables mentioned in the earlier sections, 
namely the Lorentz number and the thermopower. We begin by writing the Hubbard model Hamiltonian in the grand 
canonical ensemble 

K = ^K{r) (51) 

r 

a,fi T 

Here if represents all the neighbors of the site r, t{f]) is the hopping matrix element, and we have adopted the "midpoint 
rule" for defining densities of all non local objects such as hopping so as to get convenient expressions for the fourier 
transforms of the grand canonical ensemble Hamiltonian K = H — fiN 
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k{k) = ^e'''-^K{r) (53) 

r 

= 13 (Ep - + ^'"''^ "n^rT , (54) 

p,cr r 

with £p- = — X]^i(??) exp (— i^.p) the kinetic energy and rif? = J2a''^f,a- A more general model, with arbitrary two 
body interaction between particles, including the homogeneous electron gas (HEG) model is represented by 

k(k) = y^(e0- u) d ir + V U(q)ci ^ ir ct ^ c,- , (55) 

p,q,T,cr,cr' 

By specializing to U{q) — U and restricting the momenta to the first Brillouin Zone we recover the Hubbard model, 
whereas the electron gas is obtained by letting U{q) = 47re^/|g|^ and letting C ^ fl the box volume. In all cases, the 
electronic charge density is p{'r) = qeUp where qe is the charge of the electron (qe = — |e|). The charge current density 
is given by J(r) = Xl^o- c]f_|_i^(^c^-iif,cr) and hence its Fourier transform follows as 

p 

where the velocity vector = = ^ X^^*(^)'?exp {—iff.p). We set the lattice constant ao = 1 throughout this 

paper. 

A. The Thermoelectric operator 

Let us now work out the thcrmo electric term defined in Eq(25). The commutator can be taken and on further 
carrying out the limiting procedure we find with a = —a 



4h 



This object can be expressed completely in Fourier space as 



The more general HEG model yields the result: 



2 



qe_ 



(59) 



B. Heat Current and The Thermal Operator 6""^ 

We first derive the expression for the heat current: this follows from Eq(ll). Using Eq(54) and the fact that the 
0{U^) term vanishes ( both terms are functions of n^^s only), we find 
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2% 

u 



(60) 
(61) 



l,p.(f.a,a' 



We also need the current operator at a finite wave vector cliosen along the x axis as fc = f fc^, this may be written in 
the case of the more general HEG model as 



J'^ik) ='S^ v%isff- u) ct ir + U(q] {vf + vf, cl ^ c.- in ct ^„,Cff„' 



(62) 



We evaluate the thermal operator by a direct calculation: 



+i E m 



p,g,o- 



'^p+q,a'^P>'^P-q 



p,q,q',cr 



da' VP+9+g'^V«' ^p)^^p+q+i',a^'''P-i'^-i'- 



(63) 



Here we denote = u^(£p- — ^) and the density fluctuation = ^^ct^-^c^,;^. This expression simplifies somewhat 
in the case of the Hubbard model in real space with U{q) — > U, where it may be written as: 

d U"^ 



?7,(T 



?7 



rf,ff' ,(T 



77,77', cr 



(64) 



C. Disorder by a site potential 

We work out the various sum rules in the presence of site diagonal disorder produced by a (random) potential 
Kd = Vj' np , where we use the subscript d to denote the disorder contribution to various objects. While such a 
potential disorder does not change the electrical current operator, it does change the heat current, and it is easy to 
see that 

= ^E^-K-;+r+-|)4+.-;>V' (65) 
r,<T 
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where we defined the fouricr transform of the potential = ;^ ^-exp(ifc.r)Vf-. Prom these definitions, we find the 
disorder contribution to the thermopower sumrule 



= ^ E *(^) + ^r+^^)4+^,.Cr,a (66) 
ri,a- 



We compute the disorder contribution to the thermal conductivity sum rule as 



4 



-M X] *(^) (^f + ^r+^)4+^,aC?,'^ 

r,ff,(T 

-I E ^i^)tiv')iVx + vLf {Wr+Vr+^ + Vr+ff, + Wr+ff+f,>}cl^^^^,^^cr,a (67) 
This expression can be written in fourier space as 



l,pj' ,(7 

1,1' ,p,q,a 



(68) 



Here we need to average the final expression over the disorder potential, either perturbatively or numerically exactly. 
We save this task for a future work. 



D. Free Electron Limit and Comparison with the Boltzmann Theory 



It is easy to evaluate the various operators in the limit oi U ^ 0, and this exercise enables us to get a feel for the 
meaning of these various somewhat formal objects. We note that 



(69) 



Here rip- is the fermi function and the factor of 2 arises from spin summation. At low temperatures, we use the 
Sommerfield formula after integrating by parts, and thus obtain the leading low T behaviour; 

„2 „ ^..\/^„.x\2\ 



3 

eTT Kg 



(70) 
(71) 
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where /Oo(m) is the density of states per spin per site at the fermi level fi and the primes denote derivatives w.r.t. 
/Lt, the average is over the Fermi surface as usual. These formulas are indeed very close to what we expect from 
Boltzmann theory. Indeed replacing i/w ^ t with a relaxation time r in the leading (high frequency) terms in the 
general formulae reproduces the familiar Drude-Boltzmann results [29] for the thermal conductivity, thermopower and 
the electrical conductivity. We may form the high frequency ratios 



(72) 



It is therefore clear that the high frequency result gives the same Lorentz number as well as the thermopower that the 
Boltzmann theory gives in its simplest form. Clearly one can get more sophisticated with respect to the treatment of 
transport issues, such as a dependent relaxation time in the Boltzmann approach [28,30], provided the interesting 
Physics is in that direction. Our approach enables us to provide a more sophisticated treatment of interactions 
and disorder by computing the averages of the more complex operators given above, and this is useful for strongly 
correlated systems provided the transport relaxation issues are relatively benign. 



E. Periodic Anderson Lattice 



We present the thermoelectric operator for the important case of the periodic Anderson Lattice describing conduc- 
tion electrons c's that hybridize with a correlated set of localized levels described by f's as 

-fs^ = - X) *(^)4+77Cr - M nf<T + (eo - ti) X fl,J^^'' 

a,ff cr l<v<Nf 

+u X 4Ay+T.^(^'^)(4,jr,. + {h.c.)). (73) 

The charge current operator is exactly as it is in the case of the Hubbard model since the / levels are localized, and 
using the same ideas as before to define an inhomogeneous energy K{k) and computing the familiar commutator 
Eq(25) we find the thermo electric operator 

= E I; Hi^p - M)} 4>^- + 1 E -)^(4,.4> + flA)- (74) 

At T = the expectation of this operator should vanish (see the discussion section), whereby we obtain a formal 
expression for the chemical potential at T = from the above. 



VIII. STRONG COUPLING MODELS: 
THE HEISENBERG MODEL AND THE U = oo HUBBARD MODEL 



In this section we consider the case of strong coupling, and present the Thermal and Thermoelectric operators for 
the case of the Heisenberg model and the infinite correlation limit of the Hubbard model. 



A. Heisenberg model 



We first obtain the thermal operator for the ever popular Heisenberg model in any dimension described by the 
Hamiltonian 



H — Hp ; Hp — — JfjSf.Sf-^-ff. (75) 
f f} 

Clearly H{k) = ^ e^'^ '^Hp and from the conservation law of local energy we obtain the energy current operator 
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'r,j)i,j)2 

Using the standard commutators we find the thermal operator to be a four spin operator in general. It is expressed 

as 

= g E ('71, x ^ V2,x)JrfiJfi2'^f!3 [%,x('5f+77i X ^f+fj^).{Sf X Sf^fj^) 

-(2??2,a; + %,a;)(S'?+ifi X Sp).{Sp+fl^ X Sf^+^g) + (2r7i,x + r?3,a:)(S'?+ifi X Sp+f;^) .{Sp+fJ^ X S'f)]. (77) 

In this expression f,f+f]i,f+f]2 are necessarily distinct sites, but f+ffs can coincide with the last two, and hence 
the order of the products needs to be treated carefully. For simple cubic lattices this expression simplifies a bit, but 
in essence it involves four spin correlations. 



B. The U ^ oo Hubbard model 



In this section we consider the limit of /7 — > oo, and consider the kinetic energy only, i.e. the t part of the t — J 
model, since this is expected to dominate in transport properties, at least far enough from half filling and for f >> J. 
The addition of the J part can be done without too much difficulty, but we ignore it here for brevity. In this limit 
the fermionic commutation relations need to be modified into the Gutzwiller-Hubbard projected operators [31] 

Cf,<7 = Pg Cr,a Pg 

= Y.,a' Sp^p (78) 

The presence of the Y factor is due to strong correlations, and makes the computation nontrivial. The number 
operator np^a is unaffected by the projection. We note down the expressions for the charge current and the energy 
current at finite wave vectors by direct computation: 

K{k) = -J2 m + M%,o) e^'-^^^+ivl ct^. ^c,v 



r, 77, 77', (7 

We evaluate the thermopower operator as: 

^^x. ^ _| ^ ^ ^/ )2 ^( ^( Y^,,^{r + cl^.^~,^^,cr,a - qeH E ii^) 4+f?,.c.>. (80) 

fi,ri' ,a,a' ,r 

Note that if we set F — > 1, this expression reduces to the — > limit of the Hubbard result Eq(58), as one would 
expect. We next compute the thermal operator (summing over all spin variables) 



if,!]' ,77" ^r^(T,<T',<T" 

+ \ E + ^x) + v'. + €) m tin') t{rf') 

ff^r}' ,r)" ,r,(T 
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IX. THE TRIANGULAR LATTICE SODIUM COBALT OXIDE: HIGH TEMPERATURE EXPANSION 

FOR THERMOPOWER 



The Sodium Cobalt Oxide N axC 0O2 system is of great current interest; the composition x .68 gives a metal with 
a high thermopower ~ IQO^V/K, which is further highly magnetic field dependent [32,33]. It is notable in that the 
underlying lattice is triangular, and it has been modeled by a triangular lattice t — Jmodel with electron doping [34]. 
In this section we give a brief application of our technique to this system, which yields an interesting formula for the 
high T limit of the thermopower S* that contains significant correction to the Heikes Mott formula that is used in the 
same high T limit. A detailed study is in preparation and will be published separately, here we use the leading high 
T term to illustrate the advantage of the above formalism in tackling such a problem. We neglect contributions from 
the exchange part of the t — J model and focus on the kinetic energy which is expected to dominate the transport 
contributions. Let us compute the thermopower S* from Eqs(80,A4,29) 

= ^^^^ 

where 

^ = E + ^i)' *(^^) iYa'Ar + V) c\^^^^,^^,c,A (83) 

This is a very useful alternate formula to the Heikes formula [35,36], where the second term in Eq(82) is thrown out. 
It interpolates very usefully between the standard formulas for low temperature as well as at high temperature. The 

second term represents the "transport" contribution to the thermopower, whereas the first term is the thermodynamic 
or entropic part, which dominates at high temperature, as we shall show. In fact for S* we can actually make a 
systematic expansion in powers of unlike the dc counterpart, where it is not possible to make such an expansion. 

The computation of the different parts proceeds as follows: we show readily that (for the hole doped case) using 
translation invariance and with n as the number of particles per site at high T, 

n(T^^) = 6£g^^(ctco) ~ 3£g^/3i2n(l - n). (84) 

The structure of the term Eq(83) is most instructive. At high temperatures, for a square lattice we need to go to 
second order in (it to get a contribution with % + 7?^ 9^ 0, to the expectation of the hopping (ct^^^^, ^,Cr,cr)- For 
the triangular lattice, on the other hand, we already have a contribution at first order. For the triangular lattice, 
corresponding to each nearest neighbor, there are precisely two neighbors where the third hop is a nearest neighbor 
hop, and each of these gives the same factor (r\x + ??i)^ = 1/4 hence at high temperatures 

^ - E(^^'.<^(^l<^+V...Sg, J. (85) 



of exp(— and hence a simple estimation yields 



The spins must be the same to the leading order in fit where we generate a hopping term cl c-. , -> ^ from an expansion 

0,(T ViV i*^ 



A = -^£t^f3n{l^n){2-n) + 0{(3^). (86) 



This together with /i/fcsT = log(n/2(l - n)) + 0{pH'^) gives us the result for < n < 1 



S* = Y |log[2(l - n)/n] - (3t^ + 0{pH') \ , (87) 



and 



S* = -^ {log[2(n - l)/(2 - n)] + pt^ + 0{l3H^)] (88) 



for 1 < n < 2 using particle hole symmetry, with = — jej. Thus we can estimate the hopping parameter from the 
high temperature T dependence of the thermopower. From the temperature dependence of the data of Terasaki et 
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al( [32]) and assuming S ^ 5*onc finds that t = —llO'^K, and with this, S* ^ 120^V/K, fairly close to the observed 
value. Thus the sign and rough magnitude of the hopping are as expected from other grounds [34]. An interesting 
corollary of our analysis is that if the sign of hopping t were reversed, i.e. if t > for the electron doped case, then 
the thcrmopowcr S* would approach its high T value from, above, and since the low T behaviour is fixed to vanish 
linearly, it implies that S* (and hence presumably S), must have a maximum, unlike the case t < 0. It follows that 
a material with similar absolute value of hopping as NaxCoO^, but with the opposite sign of hopping, would have a 
greater and hence an even more exciting thermoelectric behaviour. 

Hence, the nature of the high T expansion shows that the triangular lattice is exceptional in that the transport 
corrections to the entropic part are much larger than for the square lattice; they are of 0{(ii) for the first as opposed 
to 0{(it)'^ for the latter. Also these are potentially much more sensitive to a magnetic field than the corresponding 
one for the square lattice since these 0{(5t) terms are themselves functions of the magnetic field ~ 0{B'^). In the 
degenerate limit, it seems quite possible that the field dependence of the transport part might be comparable to the 
spin entropy contribution [33]. 

The alternate formula Eq(82, 29) for the thermopower thus has the interesting property that it captures the expected 
low T fermi liquid type behaviour as well as the high T Heikes type behaviour. It should be most interesting to apply 
numerical techniques to evaluate this for all T, a calculation that seems quite possible with existing techniques, at 
least for small systems. 



X. DISCUSSION 

It is curious that the recognition of the sum rule Eq(18) for the thermal conductivity, has lagged so far behind 
Kubo's seminal paper in 1957 [10]. One of the reasons may be that most later workers followed the method given in 
Sec IV following the Kubo identity. This method does give the sum rule as shown here, but only when one dissects 
the Kubo identity carefully, keeping the possibility of superconductivity in mind. On the other hand, the many body 
linear response type method adopted in Sec II, runs into some discouragement, which is relieved only upon recognizing 
the role of the Identity I in leading to a sensible result, as elaborated in the discussion following Eq(12). 

We have shown that the sum rule for thermal conductivity involves the thermal operator 6^^ and the Seebeck 
coefficient involves which are formally evaluated in this work, for various models of current interest. We plan to 
return to a numerical evaluation of some of these in the context of strongly correlated matter in a future study. We 
make some remarks on the nature of the variables and the prospects for their evaluation. 

One notable fact is the vanishing with of the expectation of ©^^ and for a fermi gas, this is seen from 
the evaluation in Eq(71) explicitly. This leads to questions such as: What is the origin of this behaviour? Is this 
behaviour true in general? The answer is that these variables must vanish at low temperature in a fashion that is 
dictated by the specific heat. This connection is a deep one and we explain it next. 

Let us first recognize that the specific heat (at constant chemical potential) can be written as an energy-energy 
correlation function: 

r/3 



^lim / dT{K{k,-iT)K{-k)), (89) 

and the T dependence of this correlator is known in various systems in general terms. For metals and other fermi 
liquids, the specific heat at constant N, i.e. Cjv is the usual measured one and differs from this by the thermodynamic 
relation 



(as/art 



where the correction term is small for a fermi liquid at low T ~ O(T^). The thermal operator is expressed as 

= iC.vlff (91) 



~hT' ~ d^^'^^ff 



where "d" is the spatial dimension, and the above defined effective velocity Vgf / is given by 

-d d^ 



J„^rfT(jg(0,^ZT) -4^(0,0)) 

' J^dT{K{0+,-iT)k{0~,0)) 
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Here we have used the relation Dq = for generic systems to replace by the current current expectation. The 
above formulae show that the variable ^jr- is best thought of as being the specific heat times the square of a velocity, 
much as in the kinetic theory of thermal conduction where k = ^Cv^t, so that the temperature dependence comes 
predominantly from the specific heat and the effective velocity captures the dynamics, without being very sensitive 
to T. Of course it must be understood that Q^^ does not contain the scattering rate, which disappears in the sum 
rule. From a computational point of view, it seems best to use such a decomposition. At any rate, this line of 
argument gives us a qualitative understanding for the rather "magical property" of the O^^ operators listed, namely 
the vanishing of their average at low T in the true ground state. 

Our work also yields an interesting exact formula for the zero temperature chemical potential for most metallic 
many body systems. Since we expect the thermopowcr of metals to vanish in the ground state at all frequencies 
(being related to the entropy), equating the expectation of to zero gives us explicit formulae for the ground state 
chemical potential /i(0). This approach can be used in most cases, including hard core bose metallic systems. In the 
interesting case of the Hubbard model, from Eq(57), we find a formula for the ground state chemical potential as a 
ratio of expectation values of two operators 

V 

n 

where the denominator is recognized as essentially the stress tensor (r^^). In the non interacting case, this reduces 
to /i(0) = Ef-, and thus generalizes this familiar relation to the interacting case. Similar formulas result for other 
models considered in this work, using the same idea. Usually one has to resort to differentiating the ground state 
energy with respect to the number of particles, and this can be inaccurate in numerical studies. Our formula is exact 
and a consequence of the vanishing of S** at T = 0. It is interesting that standard many body text books [11,17] 
do not quote this type of expression, and hence it is interesting to verify this in some cases where exact results are 
known from other arguments. The Hubbard model on bipartite lattices is a good a check of this result, since the 
chemical potential at half filling (only) is easily found using particle hole symmetry to be U jl. One can manipulate 
the expression Eq(93) using the particle hole symmetry, and this result is easily reproduced. 

We have presented the detailed form of the thermal and thermo electric operators 9^^ and for the Hubbard 
model, and also for the strongly correlated limit /7 — > oo above. Results for other models not discussed here, can be 
readily found using the method presented. We have also presented the operators that give the disorder contribution 
to these variables, these can in principle be evaluated in perturbation theory in the disorder and interaction, and it 
should be interesting to compute these as well in low dimensions where the effect of disorder is marked. The operators 
for the lattice thermal conductivity should be interesting in the context of non linear lattices, and lend themselves 
to numerical evaluation rather easily in the classical limit, and the quantum cases also seem to be manageable with 
existing computational resources. 

One set of applications concerns the effect of strong correlations on the thermo electric power factor and the figure 
of merit in narrow band systems. These can be evaluated at high frequencies and such calculations should be useful 
guides to the role of interactions and band filling. We provide in Section IX, a simple example of this formalism by 
computing the high temperature thermopower for the triangular lattice. This important and currently popular case 
models the physics of sodium cobalt oxide [33]. Interestingly the T dependence of the thermopower leads to estimates 
of the band width that seem comparable to other estimations. 
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APPENDIX A: ELECTRICAL CONDUCTIVITY 



To complete the calculation we summarize the results for the electrical conductivity, which are already well known 

in literature [6,13,18,19,21,22]. We note that in parallel to the calculation of the thermal conductivity, the conductivity 
can be expressed in terms of the correlations of the current operator as: 



a{u)c) 



where the stress tensor 



i 1 
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e^-=*'dt'([4(t'),Jx(0)]) 
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- ct Ct. 
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dkl k,a k,a 



(A2) 
(A3) 
(A4) 



assuming a charge qe for the particles, and n{r) the particle density at r, and the charge fluctuation operator p{k) 
g'e X^exp(ifc.r)n(r). Again performing the Lehmann representation we find 



i 1 



Pn Pm 



■|(n|4|m)|2 



(A5) 



We can use partial fractions and write 



where the charge or Meissner stifl'ness is given by 



(A6) 



Dc = ^ 
C 



Pn Pm 



Kn|4|m)p 



(A7) 



The more familiar superfluid density ps{T) and plasma frequencies Wp, Wp^g for the total and superconducting conden- 
sate are defined in terms of Dc and (r^^) by 



2 _47r9,V,(T) A-kD,{T) 



h 



^2 ^ ^/^xx.^ 



(A8) 



so that Eq(A7) is just the lattice version of the well known London decomposition of electronic density into super- 
conducting and normal parts. The real part of the Lehmann rep is 



Re a{u)) = ■jTS{hco)Dc + ( ) 2^ pn\{n\Jx\m)\'^6{€m - Cn - ^<^), 



(A9) 



where 



Pn -Pr, 



{n\Ja:\m)f 



(AlO) 



using a cancellation between the terms with equal energy between the two terms. 

The lattice version of the f- sum rule for the real part of the conductivity (an even function of w) follows as 
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(All) 



We may rewrite Eq(A6) in a more compact form as 




'c '-' Jo Jo 



1 fO° fP 



(A12) 



The difference between the two stiffnesses is 



(A13) 



Thus if we know that Dc = by some independent argument, then this provides us with an alternate expression for 
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